Numerical study of the strongly screened vortex glass model in an external field 
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The vortex glass model for a disordered high- Tc superconductor in an external magnetic field is 
studied in the strong screening limit. With exact ground state (i.e. T = 0) calculations we show 
that 1) the ground state of the vortex configuration varies drastically with infinitesimal variations of 
the strength of the external field, 2) the minimum energy of global excitation loops of length scale L 
do not depend on the strength of the external field, however 3) the excitation loops themself depend 
sensibly on the field. From 2) we infer the absence of a true superconducting state at any finite 
temperature independent of the external field. 
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I. INTRODUCTION 



The gauge or vortex glass model has become a 
paradigm in studying amorphous high-Tc supierconduc- 
tors or random Josephson-junction arrays (seeo for a re- 
view). One essential feature of this model is the possi- 
ble appearance of a glassy state at low enough tempera- 
tures, without which true superconductivity (i.e. vanish- 
ing resistance) would cease to exist in these disordered 
materialsHa. 

Experimental evidence of such a vortex glass state has 
been reported for high-Tc superconductorsQ. From the 
theoretical side it is now commonly believed that in the 
absence of screening a true superconducting vortex glass 
phase occurs at low enough temperaturescTQ. If screening 
is present the original, unscreened 1/r-interaction of the 
vortex lines is exponentially shielded beyond a particular 
length scale A and the situation seems to change, in par- 
ticular in the limit in which the screening length is zero 
(i.e. where vortex lines interact only on-site) the low t^ffi- 
perature vortex glass phase seems to be destroyed0nl3tl. 

In a typical experimental situationQ the amorphous 
high-Tc superconductor is put into a homogeneous mag- 
netic field pointing, say, in the z-direction. Due to bulk 
disorder, i.e. inhomogeneities (vacancies, defects, etc.), 
in the bulk of the sample the vector potential acting on 
the superconducting phase variables attains a random 
component (most plausible this mechanisrri-J£ explained 
in the context of granular superconductors^) , however, 
still there should be a homogeneous back ground field 
superposed on the random part. 

Therefore in this paper we study the question of how 
is the latter scenario, i.e. the absence of a true supercon- 
ducting phase in the strongly screened three-dimensional 
gauge or vortex glass model influenced by the presence of 
a homogeneous external field in one particular space di- 
rection. This is done via the investigation of exact ground 
states of the vortex glass Hamiltonian and its low en- 
ergy excitation. First we analyze the sensibility of the 
minimum energy configuration with respect to the addi- 



tion of a homogeneous external field, then we study the 
low energy excitations of length scale L in the spirit of 
the usual domain wall renormalization group (DWRG) 
calculationsQOH. 

The lattice model describing the phase fluctua- 
tions (described by phase variable (j)i G [0, 2tt[) in a 
strongly disordered superconductor close to a normal- 
to-suDcpconductor phase transition is the gauge glass 
modeU'El 
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where the first sum runs over all nearest-neighbor pairs 
{ij) on a L X L X L simple cubic lattice and the sec- 
ond over every elementary plaquette □, respectively, and 
the fluctuating vector potentials and A the screening 

are the quenched vector po- 
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tentials consisting of a random component Af- £ [0, 27r[ 
and a homogeneous component A^°™ modeling an ex- 
ternal magnetic field in the z-direction. The parameter 
A is the bare screening length. A similar Hamiltonian 
occurs for ceramic (granular) supempnductors including 
the self-inductance of vortex loopsll3"li3. For simplicity 
we set J = 1 . n n 

After standard manipulationalj'EZl this Hamiltonian 
can be brought into the vortex representationlS 
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The Ji, representing the vortex density of the phase fields 
(j)i on bond i, are three-component integer variables run- 
ning from — CX3 to oo, live on the links of the dual lattice 
and satisfy the divergence constraint (V- J)^ = on every 
site i. The function G(r) is a lattice Greens function be- 
having asymptotically like G(r) ~ r~^ exp(— r/A). The 
bi are magnetic fields which are constructed from the 
quenched vector potentials Aij by a lattice curl, in the 
present case with the homogeneous external field we have: 
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Obviously the random part fulfills the divergence free 
condition (V • b)^ = 0. We specify the boundary condi- 
tions for the vortex glass Haniiltonian to be periodic in 
all space directions (corresponding to fluctuating bound- 
ary conditions in tiie jihase variables of the original gauge 
glass Hamiltonianll3^tj. Now we choose B"'^* = Bez i.e. 
the external field points in the z-direction (i.e. along e^ 
the unit vector in the z-direction) and is also divergence- 
less due to the periodic boundary conditions. 

We just remark that in the pure case (A'''"'^ = 0) the 
field strength B simply plays the role of the usual fill- 
ing factor / counting the number of flux units per pla- 
quette gpiipe rise to the uniformly frustrated XY-model 
(see e.g.E3'El and references therein) in the unscreened 
case (A — oo). Here, due to the long range interaction 
G{r) « l/r, the ground state is indeed non-trivial for 
irrational filling factors. In the continuum limit the flux 
lines would actually form a hexagonal lattice, the well 
known Abrikosov flux line lattice. 

For the disordered case one has an interesting interplay 
between two sorts of frustration: one is also present in 
the pure case and coming from the external field and the 
other comes from the quenched disorder. To our knowl- 
edge this problem has not been investigated systemati- 
cally so far. Here, as a first step, we confine ourselves to 
the strongly screened casc_(A -^ 0), for which the vortex 
Hamiltonian (0) simplifiesO to 



H, 
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The problem of finding the ground state, i.e. the min- 
imum energy configuration of this Hamiltonian, is ac- 
tually a minimum-cost-flow- problem minjj.j. ^^ Cj(Ji) 
subject to the constraint (V • J); = 0, where 
Ci(Ji) •— (Jj ~ bi)^/2 are the so called (convex) cost 
functions. This problem can be solved exactly in polyno- 
mial time via, combinatorial optimization techniques, as 
described inBO. 



II. FIELD CHAOS 



First we investigate how the ground state of the vor- 
tex glass changes with the external field. Obviously, 
with increasing field strength B more and more flux 
will flow in the z-direction. In Fig. 1 we depict var- 
ious quantities reflecting this observation: a) the total 
flux L~^Yl,ii\-^i\ + \Ji\ + I'^fl): which increases, after a 
crossover region, linearly with B; b) the fraction of flow 
variables in the x- and y-direction that are zero/nonzero; 
then the fraction of flow variable in the z-direction which 
are c) larger than, d) smaller than, e) equal to zero. 
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FIG. 1. Plot of the count of components per site 

vs. external field for 1000 samples. a) The total flux 

^-3 v^ r\ Tx\ , \ jv\ , \ rz|^. u^ r-3" 

L- 



'EM'l + \Jf\ + \J^\)^ b) L-^E.'5./r,o (upper line), 
'Y.^^(+Ji) and L-^Y.r^i-Jf) (both lay on the lower 
line); the j/-values are the same; c) L""^ ^. S(-|-Jf ); d) 
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For B — Q the vortex variables Ji are homoge- 
neously distributed; L~'^Yl,i^Ji,Oi L~^J2i^{~^Ji) ^^^ 
L~^'Y^^0{—Ji) have the same values for each direction, 
where the last two values are equal due to the sym- 
metry of direction. An increasing field B > Q only ef- 
fects the components along the z-direction (dashed and 
dotted lines) and let the components perpendicular to 
it unchanged (solid lines). There is a critical external 
field Be, above which all components of the flux lines Jf 
parallel to the field have a non-vanishing negative am- 
plitude, what is attended by a decreasing of the zero- 
components, i.e. for B -^ Be the #{Jf < 0} increases 
(Fig. 0d) and the #{Jf = 0} decreases (Fig. |c). Be- 
fore all zero-components vanish {Be = 5/2), the positive 
components have disappeared (Be = 3/2). The criti- 
cal fields can be determined by remembering that the 
random field b*""" is in [—2, -1-2] and the fiux is an inte- 
ger. Thus the maximal discrepancy between the optimal 
flux and field is 1/2 per site. The optimal Jf depends 
on the applied field by a step function and in the worst 
case the random field is b™" = ±2, so that we have 
{J°P*} = {J, integer] min |J, - (B^^* ± 2.0)| }. One 
obtains: Jf < Vi if B > B^, = 5/2 and Jf < Vi if 
B > Bc2 — 3/2, in accordance with Fig. 1. 

Next we study the sensibility of the ground state (or 
optimal flow configuration) with respect to small changes 
in the external field B. To this end we compare the 
ground state configurations of samples with the same 
quenched disorder and slightly different external field B. 
Denoting with Ji the zero field (J5 = 0) ground state 
and with Ji(-B) the ground state of the same sample in 
non-vanishing external field B we define the Hamming 
distance of the two configuration J^ and Ji (B) by 
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FIG. 2. Scaling plot of the Hamming distance Db{L) vs. 
LB^/'^ for L < 32: 5000 samples for L < 16, 2000 for L = 24 
and 500 for L = 32. The chosen values for B are B = 0.0001, 
0.0010, 0.0100 and 0.1000. The best data collapse is achieved 
by a chaos exponent (^ = 3.8 ± 0.2. The error bars are less 
than the size of the symbols and thus omitted. 



DB(i)=^(J,(B)-J, 



(5) 



so that a small value of Db (L) means a strong correlation 
of tha.ground states. 

Ir£3 it has been found that an infinitesimal random 
perturbation of the vector potential AJ^"'* leads to a 
chaotic rearrangement of the ground state configuratiaa. 
There it was demonstrated that, like in spin glassesc^O, 
beyond a particular length scale, the so-called overlap 
length, the two ground states (perturbed and unper- 
turbed) decorrelate. For the non-random magnetic field 
perturbation we study here we take over this concept and 
demonstrate the existence of an overlap length I* scaling 
with the strength of the external field B like I* oc B^^'' , 
where ( is the chaos exponent. For this length scale to 
exist the finite size scaling form Db{L) = d{L/l*) ~ 
d{L/B^/'') should hold, which is indeed satisfied, as is 
shown in Fig. 2. Wc obtain a relatively large chaos ex- 
ponent C = 3.8 ± 0.2. Remarkably this exponent coin- 
cides (within the error bars) with the chaos exponent for 
a random perturhation which has been reported to be 

^rand ^ 39 ^ q j^ 



III. DEFECT ENERGY (DWRG) 

In this section we study the scaling behavior of low en- 
ergy excitations AE{L) of length scale L (to be defined 
below) in the presence of an external field, which provides 
the essential evidence about the stability of the ground 
state with respect to thermal fiuctuations. If AE{L) de- 
creases with increasing length L it implies that it costs 
less energy to turn over larger domains thus indicating 



the absence of a true ordered (glass) state at any non- 
vanishing temperature. Usually one studies such exci- 
tation of length scale L by manipulating the boundary 
conditions (b.c.) foij-|the phase variables of the original 
Hamiltonian dl^), secBu. One induces a so called domain 
wall of length scale L into the system by changing the 
b.c. of a particular sample from periodic to anti-periodic 
(or vice versa) in one space direction and measures the 
energy of such an excitation by comparing the energy of 
these two ground state configurations. This is the com- 
mon procedure for a DWRG analysis, which, however, 
bears some technical complicationsB and some concep- 
tual ambiguitiesBt] in it. 

Here we follow the basic idea of DWRG, we will, how- 
ever, avoid the complications and the ambiguities that 
appear by manipulating the b.c. and try to induce the 
low energy excitation in a-different way, as it has first 
been done by one of us ii£3 for the zero-field case. 

First we clarify what a low energy excitation of length 
scale L is: in the model under consideration here it is 
certainly a global vortex loop encircling the 3d torus (i.e. 
the L X L X L lattice with periodic b.c.) once (or sev- 
eral times) with minimum energy cost. For the pure case 
the global minimum energy loop is simply a straight line 
that costs energy AE{L) = JL, which is exactly also 
what one would expect for a domain wall of length L 
in a three-dimensional XY-model, and which is also ob- 
tained from the energy difference between ground states 
with periodic and anti-periodic b.c. 

Next we have to clarify how we induce the above men- 
tioned global vortex loop, if not by manipulating the 
b.c. The solution is to manipulate instead the costs 
for flow in one particular space direction. Suppose we 
have the exact ground state of one particular sample, 
which specifies also the current cost for increasing the 
the flow variables J^''^'^ with respect to J° by one unit, 
e.g.: Acf = c^{J°'' + 1) - c,{Jf'') = J°^ - &f + 1/2. 
If we decrease smoothly these variables Acf and apply 
our min-cost-flow algorithm to this modifled problem at 
some point a configuration 3\ that is the original ground 
state plus a global loop in the x-direction appear as the 
new optimal fiow configuration for the modified problem. 
This extra loop, which can be easily identified by com- 
paring the new optimum with the original ground state, 
is the low energy excitation we are looking for. Its en- 
ergy AE{L) is simply the difference of this state (ground 
state with loop), H{J^) minus the energy of the ground 
state _ff(J*'). Note that this is always positive, since it is 
definitely an excitation (in contrast to the usual DWRG 
procedure where the b.c. is modified). 

Four remarks are in order: 1) small, topologically sim- 
ply connected loops are not generated by this procedure, 
since all what can be gained in energy is lost again on 
the return. 2) In the pure case this procedure would not 
work, since at some point spontaneously all links in the 
z direction would increase their flow value by one. It is 
only for the disordered case with a continuous distribu- 
tion for the random variables b^ that a unique loop can 



be expected. 3) Sometimes (in ca. 5% of the samples) 
the global flux changes discontinuously by more than one 
unit, a typical example for such an elementary excitation 
loop is shown in Fig. ^, however we define these still 
to be elementary excitations of length scale L. 4) In the 
presence of a homogeneous external field one has to dis- 
criminate between different excitation loops: those par- 
allel and those perpendicular to the external field need 
not to have the same energy (however, it turns out that 
the disorder averaged defect energy is identical in all di- 
rections, see below). 

Schematically the numerical procedure is the following: 

1. Calculate the exact ground state configuration {J°} 
of the vortex Hamiltonian(0); 

2. Determine the resulting global flux along, say, the 

X-axis /:r = rE^-^^; 

3. Study a minimum-cost-flow problem in the actual 
cost for increasing the flow in the x-direction Acf = 
Q(J°^-hl)-c,(J^^)-6f-|-I/2 is smmothly modified 
letting the cost of a topologically simple connected 
loop unchanged and only affecting gloabal loops. 

4. Reduce the Acf until the optimal flow configura- 
tion {J^} for this min-cost-flow problem has the 
global flux {fx + 1), corresponding to the so called 
elementary low energy excitation on the length 
scale L; 

5. Finally, the defect energy is /S.E — _ff({J^}) - 
i7({J"}). 



It turns out that the computation time grows linearly 
with B and it took about 150 seconds per samples for 
i = 24 with B = 1 on a SUN workstation (167MIIz). 

As for the zero-field casetEl one expects for the disor- 
dere averaged excitation energy (or defect energy) 



[AE{B,L)l^^L'' 



(6) 



where B is fixed, [• • -Jav denotes the disorder average and 
9 is the stiffness exponent and its sign determines whether 
there is a finite temperature phase transition or not, as 
explained above, li 9 < 0, i.e. the transition to a true su- 
perconducting vortex state appears only at zero temper- 
ature T = 0, scaling predicts that the thermal correlation 
lengtliX, diverges at T = as ^ oc T'" with v = l/\9\, 

c.f.llia'il. 



In the non-random case A; 



rand 



an elementary exci- 



ation of the kind we described, i.e. a global loop of length 
L around the torus costs an energy of AE ^ i+^ with- 
out any dependence of the external field as well parallel 
as well as perpendicular to the external field (i.e. the z- 
direction). Note, that this is the situation for of the pure 
classical 3d XY-model with strong screening. 

For a single configuration we find that a change of the 
external magnetic field B drastically affects the defect 
energy AE (Fig. m. AE is a piecewise linear function 
that behaves in particular intervals [Ba, Bb] as 



dAE 
dB 



2Ln^ 



(7) 
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FIG. 3. The minimum energy global exciation loop perpendicular to the external field in the z-direction is shown for one 
particular sample (L = 24) and three different field strengths B (note the periodic b.c. in all space directions), a) (left) 
B £ [0.0065, 0.0069] is in a range, where the defect energy AE varies linear with respect to the field (see inset of Fig. m. Note 
that the loop has also winding number n^ = 1 in the direction parallel to the external field. Hence dAE/dB = 2L. b) (middle) 
The same sample as in (a) with B G [0.0070, 0.0075]. In this interval the defect energy is constant, no loop along the direction 
of the applied field occurs, c) (right) The same sample as in (a,b) with B € [0.0076,0.0081]. The system is very sensible to 
the variation of applied field AS. Even for a small change by AS = 0.0001 the form of the excitation loop changes drastically 
(compare with a and b). 
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FIG. 4. Defect energy AE vs. applied magnetic field B 
for one particular disorder configuration. The field varies be- 
tween and 0.05 times one flux unit and the system size is 
L — 24. The inset enlarges the region that is studied in Fig. 
P^ in more detail. 

which can be understood as follows: the external field 
varies continuously and the interger valued fiow changes 
only in discrete steps, thus the minimum energy ex- 
ciation loop may not change in a whole interval say 
[Ba,Bt,]. In this interval AE changes linearly with B 
since H{3^) — .ff(J°) is simply proportianl to length of 
the exciation loop in the z-direction, which is n^ • L, with 
n^ the winding number of the loop in the z-direction 
(n^e {...,-2,-1,0,+!, +2,...}). 

Furthermore not only the ground state itself is ex- 
tremely sensible to the small variations of the external 
field strength (as we have seen further above), but also 
the excitation loops themselves change their form dra- 
matically, as it is exemplified in Fig. g. Only small parts 
of the loop seem to persist over a significant range of the 
field strength, see for instance the in the vicinity of the 
plane z = 20 in Fig. 13. 

Now we study the behavior of the disorder averaged 
energy of excitation loops perpendicular to the applied 
magnetic field along, say, the z-direction (full diamonds 
in Fig. 0). Note that it is only necessary to study the 
situation B G [0, 1], since all physical properties of the 
vortex glass Hamiltonian (m are periodic in the strength 
of the external field B, i.e. the filling factor. As can eas- 
ily be seen in Fig. || the defect energy [AE{B, L)]av is 
independent of the value of B. For any fixed value of B 
the finite size scalingpMlation (pi) is confirmed and gives 
9^ -0.95 ±0.04, c.iM. 

We want to note that this behavior of excitations per- 
pendicular to the applied field depends neither on the 
length of the system in the z-direction nor on the topol- 
ogy in this direction: we also studied the situation in 
which the vortex Hamiltonian (0) lives on a lattice with 
free instead of periodic b.c. in the z-direction. In this 
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FIG. 5. Linear plot of the external magnetic field vs. the 
defect energy [AE]av for L = 4, 6, 8, 12, 16, 24 (top to button). 
The elementary excitation parallel □ and perpendicular full 
diamond to the external magnetic field B^"^*. For each plot- 
ted point the number of samples varied between 20000 for the 
smallest sizes to 1000 for the largest sizes. The error bars of 
the excitations are less than the size of the symbols and thus 
omitted. 

case the external field has an appropriate source and 
sink outside the system. We find here the same result as 
before: the disorder averaged defect energy is indepen- 
dent of the strength of the external field B. 

Finally, we also observe a field-independent domain 
wall energy for elementary excitation loops (open boxes 
in Fig. Ig) parallel to the external field. Thus we conclude 
that the disorder averaged defect energy is also indepen- 
dent of the direction of the homogeneous external field. 



IV. SUMMARY 

We have studied the three-dimensional vortex glass 
model in the strong screening limit in the presence 
of a homogeneous external field. The ground state 
is extremely sensible to small external field variations. 
Ground state configurations at different field values B 
and B + AB decorrelate beyond the overlap correlation 
length I* ^ AB^'^, where C is the so called chaos expo- 
nent which we estimate to be C = 3.8 ± 0.2. This value 
agrees within the error-bars with the chaos exponent for 
random perturbations, of the quenched disorder that have 
been investigated inE3. 

For individual disorder configurations the change of 
the defect energy AEl{B) with respect to the appHed 
field B is piecewise linear, analytic and accompanied by 
a drastic deformation of the minimum energy global ex- 
citation loop. On the other hand the disorder averaged 
value of the defect energy [A£'i(i3)]av is independent of 
the strength of the external field B. Moreover it turned 



out that also the excitation loops parallel as well as per- 
pendicular to the external field yield the same disorder 
averaged value [A£^i(i3)]av Thus the scaling behavior 
of the defect energy is independent [Of B, i.e. iidentical 
to the case B = already studied inll3. Therefore, as in 
the B = case, we infer the absence of a true supercon- 
ducting low temperature phase in the strongly screened 
vortex glass model in an external field. The stiffness ex- 
ponent is again, as in the B = case, 9 = —0.95 ± 0.04, 
giving an estimate i/ = l/\0\ = 1.05±0.05 for the thermal 
correlation length exponent. 

Concluding we would like to note that it would be in- 
teresting to perform the same analysis for non-vanishing 
screening length and for the unscreened case, where due 
to the long range repulsion of the vortex lines important 
new physics might appear. In particular it is an open 
whether a homogeneous external field has a significant ef- 
fect on the existence of the low temperature vortex glass 
phase in the unscreened case. 
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